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ABSTRACT 
The  surface  of  a  large,  anechoic  tank  was  excited  to  produce  several 
different  anisotropic  sea  and  swell  surfaces  characterized  by  "typical" 

but  distinctive  frequency  spectra  and  wave  height  probability  densities. 

2  2   2 
Surface  roughness  for  specular  scatter,  g  (=4k  a  sin  e),  was  varied 

from  zero  to  four  by  selection  of  underwater  sound  propagation  constant, 
k,  RMS  wave  height,  a,  and  grazing  angle,  e.  Concurrent  measurement 
was  made  of  20  second  averages  of  squared  surface  wave  height,  and 
coherent  specularly-scattered  sound  intensity.  For  g  <_  1 . ,  coherent 
specular  scatter  from  all  surfaces,  including  pure  swell  (with  a 
strongly  non-Gaussian  wave  height  probability  density)  was  found  to  vary 
as  e"9,  as  previously  predicted  for  Gaussian  surfaces  only.  For  certain 
combinations  of  e,  a,  and  sound  wavelength,  A,  coherent  components 
significantly  larger  than  predicted  were  noted  in  both  cross-wind 
and  down-wind  scattering  for  g  >  l.j  this  is  believed  to  be  a  previously- 
hypothesized  interference  phenomenon.  Cross-sea  scattering  consistently 
produced  a  slightly  stronger  coherent  component  of  intensity  than 
scattering  in  the  direction  of  the  sea. 
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I.   INTRODUCTION 

The  last  ten  years  have  witnessed  the  development  and  deployment  of 
low-frequency,  high-power,  surface  ship  sonars.  Intelligent  use  of 
these  long-range  sensors,  and  the  ability  to  predict  their  performance 
under  a  given  set  of  conditions,  require  quantitative  understanding  of 
the  sound  propagation  involved.  Both  surface  and  seabed  reflections  of 
the  beam  can  greatly  affect  sonar  performance;  quantitative  analysis  of 
"bottom-bounce"  phenomena  generally  has  not  proved  fruitful  because  of 
the  dearth  of  knowledge  concerning  the  geography,  topography,  composi- 
tion, and  acoustic  properties  of  the  ocean  bottom.  Analysis  of  surface 
scattering,  however,  is  more  tractable  because  of  the  comparative  ease 
with  which  the  properties  of  the  reflecting  boundary  may  be  studied. 
In  addition,  surface  scattering  experiments  are  easily  scaled  to  labora- 
tory dimensions. 

When  a  beam  of  underwater  sound  is  incident  upon  the  water's  surface 
it  is  scattered,  the  nature  of  this  scattering  being  dependent  upon  the 
characteristics  of  the  interface.  If  the  surface  is  smooth,  the  scat- 
tering will  be  purely  specular.  Conversely,  if  the  surface  is  rough, 
sound  will  be  scattered  in  all  directions,  with  signal  strength  being  a 
function  of  the  angle  of  observation.  Rough-surface  reflection  produces 
a  scattered  sound  field  which  contains  both  coherent  and  incoherent  pres- 
sure components  relative  to  the  incident  pressure.  The  behavior  of  the 
coherent,  scattered  component  is  of  prime  importance  because  of  the  in- 
creasingly more  common  use  of  correlation  detectors  in  sonar  systems. 

Scheible  and  Fowler  (1968)  studied  SDecular  scatter  of  underwater 
sound  from  a  wind  driven  model  sea  at  the  Naval  Postgraduate  School; 


their  investigation  included  both  the  amplitude  distribution  and  the 
frequency  spectrum  of  the  scattered  sound. 

The  work  reported  herein  is  essentially  an  extension  and  refinement 
of  the  work  done  by  Scheible.  An  analogue  correlator  (which  Scheible 
did  not  have)  was  used  to  isolate  the  coherent  component  of  the  scat- 
tered sound  pressure o  As  will  be  explained  in  a  later  section,  this 
isolation  permitted  the  use  of  much  "rougher"  water  surfaces  for  scat- 
tering. The  types  of  scattering  surfaces  used  were  increased  to  include 
pure  swell,  sea-swell  combination,  and  wind  driven  seas  of  three  differ- 
ent heights.  In  addition,  scattering  both  parallel  and  perpendicular  to 
the  sea  direction  was  studied. 

For  all  of  these  "seas"  the  objective  has  been  to  use  the  model  ex- 
periment to  study  cases  that  are  beyond  the  capability  of  theorists,  and 
beyond  the  practicality  of  ocean-acoustic  experiments. 
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II.  THEORETICAL  BEHAVIOR  OF  THE  SCATTERED  SOUND 

A.  THE  SPECULARLY  SCATTERED  INTENSITY 

The  following  treatment  briefly  summarizes  the  statistical  approach 
of  Medwin  [Ref.  1],  which  in  turn  follows  the  development  of  Beckmann 
[Ref.  2].  The  measure  of  the  roughness  of  the  water  surface  from  which 
the  sound  is  scattered  is  defined  by  the  roughness  parameter  g. 


g     =     [2tt  -     (cos  <j>1     +     cos  <j>2)] 


where:     a  =  Root-mean-squared  (RMS)  Wave  Height 

A  =  Wavelength  of  sound  in  water 

eh  =  Angle  of  incidence  of  sound  upon  the  surface 
(measured  with  normal  to  surface) 

<j>«  =  Angle  of  reflection  of  the  sound 
For  specular  scatter,  $,  *  $«  =  ♦  and  g  becomes 

g  =  k2a2(2  cos  *)2  where:  k  =  |i  *  *| 

The  general  scattering  problem  is  then  attacked  by  applying  the 
Helmholtz  Integral  to  a  volume  within  the  water,  part  of  the  surface  of 
which  is  the  ensonified  surface  area.  If  the  Kirchhoff  assumption  (the 
use  of  the  plane  wave  reflection  coefficient  at  each  point  on  the  sur- 
face) is  satisfied,  an  expression  for  the  time-averaged,  reflected 
intensity  relative  to  mirror  surface  intensity  results,  which,  for  the 
specular  direction  reduces  to: 

irouahX   .   e-g  +  f(AiLi+,0iX) 

'  mirror/ 
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where   L  =  Correlation  distance  on  the  rough  surface 
A  =  Ensonified  surface  area 
I  =  p*p  where  p  =  acoustic  pressure 

The  above  expression  was  derived  for  a  surface  correlation  function  of 
Gaussian  form,  and  for  a  Gaussian  distribution  of  surface  wave  heights. 

The  first  term,  describing  the  coherent  specular  component,  is  iden- 
tical to  that  obtained  by  several  authors.  The  second  term  represents 
scattered  sound  which  is  incoherent  with  the  transmitted  signal;  since 
it  was  not  possible  to  experimentally  isolate  the  wery   small  incoherent 
component,  the  exact  form  of  the  term  is  not  given  here.  However,  the 
following  theoretical  conclusions  are  considered  worthy  of  note: 

1.  The  magnitude  of  the  incoherent  term  is  a  function  of  ensonified 
area  and  correlation  distance,  while  the  coherent  term,  e"g,  is 
independent  of  these  parameters . 

2.  In  the  mirror  surface  limit  (g-*0)  the  incoherent  term  vanishes, 
and  the  coherent  component  approaches  unity. 

Scheible  [Ref.  3]  used  this  limiting  behavior  to  suppress  the  incoherent 
term;  he  limited  g  to  0  <  g  <  0.3,  stated  that  the  first  term  was  much 
greater  than  the  second,  and  dropped  the  incoherent  term. 

Attacking  the  problem  through  the  surface  wave  spectrum,  Parkins 
[Ref.  4]  also  predicts  a  coherent  term  of  the  form  e~^0  His  incoherent 
term  arises  from  selective  interactions  between  appropriate  surface  wave 
and  sound  wavelengths.  In  the  purely  specular  direction,  zero  incoherent 
component  is  predicted,  since  his  treatment  considers  only  the  central  ray 
and  neglects  scattering  contributions  from  off -axis  parts  of  the  ensonified 
area. 
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B.  SURFACE  SHADOWING  AT  NEAR  GRAZING  INCIDENCE 

As  the  grazing  angle  (complement  of  angle  of  incidence,  <j>)  approaches 
zero,  shadowing  of  certain  areas  of  the  scattering  surface  will  result. 
A  shadowing  correction,  derived  from  ray  acoustics  and  a  statistical 
representation  of  the  surface,  is  presented  by  Wagner  [Ref.  5],  who 
states  that  the  correct  use  of  the  shadowing  function,  S(e),  for  the 
coherent  scattered  component  is: 


^coherent)   =   $2(e)  e-g 


mir 


where:     S(e)  =  (1  -  e"4B)  erf(^) 


B   =  J.  [(3.V)-1/2  e"(9v2)/8  ♦  M'^2   e-v'  -  erfc  v] 


tan  0 


2.(o/l) 

a      =  RMS  wave  height 

L   =  Surface  correlation  distance 

e   =  Grazing  angle 

Diffraction  effects  are  neglected  by  the  ray  treatment;  it  is  therefore 
expected  that  the  accuracy  of  Wagner's  correction  will  increase  with 
increasing  -——  . 
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III.   EXPERIMENTAL  PROCEDURE 

A.  THE  ANECHOIC  TANK 

The  scattering  experiments  were  carried  out  in  the  north  anechoic 
tank,  filled  with  fresh  water,  and  of  dimensions  24'  long,  6'  wide,  and 
8'  deep,  A  wall  and  bottom  lining  of  acoustically  absorbant  rubber 
material  renders  the  tank  essentially  anechoic  at  frequencies  above 
20  kHz.  A  swell  maker  and  wind-wave  generator  were  located  at  the  east 
end  of  the  tank  as  shown  in  figure  1.  Clear  plastic  covers,  suspended 
6"  above  the  water,  confined  the  wind  generator's  output  to  the  surface, 
producing  a  wind-driven  sea.  Five  separately  controllable,  centrifugal 
fans  feeding  a  common  plenum  provided  the  flexibility  needed  to  produce 
RMS  wave  heights  from  0.1  to  0.5  cm.  Swell  waves  were  produced  by  sinu- 
soidally  oscillating  a  "V"  shaped  wooden  wedge  in  the  vertical  plane  by 
means  of  a  motor-driven  eccentric,  as  shown  in  figure  2. 

The  sea,  swell,  or  sea-swell  combination  propagated  down  the  tank, 
past  the  sound  scattering  area,  and  was  largely  dissipated  by  a  sloping 
"beach"  consisting  of  loosely  packed  aluminum  shavings  within  a  plastic 
screen  enclosure.  In  all  cases  the  ensonified  area  was  centered  on  the 
wave  staff,  which  was  located  at  a  fetch  distance  of  3.32  meters  meas- 
ured from  the  plenum  discharge. 

The  surface  area  ensonified  by  the  active  transducer  can  be  computed 
geometrically  by  assuming  no  ensonification  outside  the  -3dB  angles.  The 
resulting  surface  is  an  ellipse,  whose  area,  A,  and  largest  linear  dimen- 
sion, d,  are  computed  as  shown  in  figure  3,  Attention  is  invited  to  the 
distinction  between  grazing  angle,  e,  and  its  complement,  angle  of  inci- 
dence, <|>.  The  grazing  angle,  e,  was  determined  by  measuring  transducer 

separation  and  transducer  depth, 
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To  scatter  sound  both  parallel  and  perpendicular  to  the  sea  direc- 
tion, three  separate  transducer  arrangements  were  used  as  shown  in 
figure  4.  The  transducers,  supported  by  the  mounting  racks  described 
in  Ref.  1,  were  always  positioned  with  their  axes  in  the  specular  direc- 
tion. For  frequencies  up  to  and  including  120  kHz,  two  EDO  model  327 
piezoelectric  transducers  were  employed;  the  face  of  each  consisted  of 
two  concentric,  electrically  separated  driving  sections.  At  a  given 
frequency,  two  beam  patterns  were  available:  a  broader  beam  obtained 
by  driving  the  inner  section  only,  and  a  narrower  beam  resulting  from 
driving  both  sections  in  electrical  parallel. 

B.  THE  SEA  SURFACE 

Various  surfaces  consisting  of  "sea"  alone,  "swell"  alone,  and 
"sea-swell"  combination  were  used  in  the  scattering  study.  To  determine 
the  nature  of  each  scattering  surface  the  instantaneous  wave  height, 
z(t),  at  the  center  of  the  ensonified  area  was  continuously  monitored  by 
an  omnidirectional  wave  staff;  the  resulting  wave  record  was  tape  recorded. 
Upon  completion  of  a  particular  scattering  study,  the  corresponding  wave 
data  were  replayed  and  processed  to  determine  the  RMS  wave  height  which 
existed  at  the  time  of  the  sound  scattering, 

1.  The  Wave  Staff 

The  wave  staff  and  its  associated  electrical  circuitry  was  de- 
signed, fabricated,  and  calibrated  locally.  The  staff,  shown  in  figure  5, 
is  of  the  resistance  type  and  operates  as  one  leg  of  a  balanced  resistance 
bridge,  as  shown  in  figure  6. 

Over  a  limited  range  of  immersion  the  resistance  between  the  staff 
and  ground  rod  is  inversely  proportional  to  the  depth  of  immersion,  and  it 
is  within  this  linear  range  that  the  staff  is  made  to  operate,  as  shown  in 

figure  9A. 
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Preliminary  experimentation  revealed  that  the  position  of  the 
ground  rod  relative  to  the  staff  affected  the  static  calibration  of  the 
instrument.  To  eliminate  this  variable,  and  to  insure  its  insensitivity 
to  passing  waves,  the  ground  rod  was  mounted  on  the  same  support  bracket 
as  the  staff  and  was  insulated  by  a  plastic  sleeve  where  it  pierced  the 
water's  surface  (figure  5). 

The  wave  staff  assembly  formed  a  time-varying  resistance  whose  mean 
value  was  about  2500  ohms;  the  other  three  legs  of  the  nearly-balanced  re- 
sistance bridge,  shown  in  figure  6,  were  2500  ohm  resistors.  This  bridge 
was  fed  with  AC  voltage  at  25  kHz  to  prevent  the  formation  of  electrolysis 
bubbles  on  the  staff,  which  would  have  caused  a  change  in  resistance.  The 
bridge  output,  a  25  kHz  voltage  with  a  time-varying  envelope,  was  amplified, 
demodulated,  and  reduced  to  zero  mean  by  a  high-pass  RC  filter  with  half- 
power  point  of  0.05  Hz.  The  output  of  the  circuit  is  called  the  wave  record, 
z(t),  a  time-varying  voltage  proportional  to  the  instantaneous  wave  height. 

a.  Wave  Staff  Calibration 

The  data  of  figure  9A  (from  which  was  computed  the  static  probe 
calibration)  were  obtained  by  measuring  changes  in  probe  circuit  output 
voltage  as  the  water  level  in  the  tank  was  changed  by  accurately  measured 
increments;  it  was  necessary  to  take  these  readings  with  the  DC  blocking 
capacitor  of  figure  6  by-passed. 

In  addition  to  this  static  deflection  calibration,  the  wave 
staff  (and  its  associated  circuitry)  was  calibrated  for  surface  wave  fre- 
quencies up  to  30.0  Hz.  At  each  frequency  checked,  an  independent 
determination  of  the  calibration  was  made.  Each  was  then  normalized 
(figure  9B)  by  dividing  by  the  previously  obtained  static  calibration. 
When  the  resulting  data  were  plotted  as  shown  in  figure  9B,  the  frequency 
response  curve  for  the  probe  and  its  circuitry  results.  In  the  frequency 
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range  0.5  <f  <  8.0  Hz,  motion  of  the  water  relative  to  the  staff  was 
obtained  by  vertically  oscillating  the  staff  assembly  in  a  beaker  of 
"tank  water",  as  shown  in  figure  7.  The  frequency  and  amplitude  of  this 
displacement  were  both  adjustable  and  accurately  measurable. 

In  the  frequency  range  8.0  <_  f  <_  30.0  Hz,  change  in  probe 
immersion  was  accomplished  by  generating  capillary  waves  on  the  surface 
of  a  small  tank.  Referrring  to  figure  8,  the  generating  mechanism  was 
a  lightweight  fiber  cone,  mounted  point  down,  partially  beneath  the  sur- 
face as  shown.  The  cone  was  driven  sinusoidally  in  the  vertical  plane 
by  an  electric  vibrator  (a  permanent  magnet,  moving  coil  device)  mounted 
rigidly  above  the  tank.  The  frequency  and  amplitude  of  the  cone's 
motion  could  be  varied  by  changing  the  corresponding  properties  of  the 
driying  voltage.  To  prevent  capillary  reflection  from  the  walls  of  the 
rectangular  tank,  a  sloping  "beach"  of  blotting  paper  was  installed 
along  all  four  sides.  Capillary  height  was  measured  both  by  the  wave 
staff  and  by  a  micrometer  assembly,  mounted  rigidly  to  the  frame  of  the 
tank,  by  means  of  which  a  thin,  pointed  rod  could  be  lowered  to  touch  the 
water's  surface.  Measurement  of  maximum  wave  height  by  the  micrometer 
assembly  was  done  as  follows: 

1.  With  the  surface  in  a  mirror  condition,  the  rod  was 
slowly  and  carefully  lowered  until  its  point  just  touched  the  surface. 
Reading  the  micrometer  dial  gave  a  surface  height  reading  accurate  to 
1.02  mm;  to  reduce  this  uncertainty,  five  consecutive  readings  were  taken 
and  averaged. 

2,  Capillary  waves  of  the  desired  frequency  were  then  gener- 
ated, and  the  above  measurement  repeated.  The  point  of  first  contact  in 
this  case  was  the  height  of  the  crests  of  the  waves. 
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3.  Capillary  height  (crest-to-trough)  was  then  computed  as 
being  twice  the  difference  in  the  two  measurements.  This  step  assumes 
that  crest-to-trough  height  is  equal  to  twice  the  value  of  mirror-to- 
crest  height. 

The  crest«to-trough  voltage  variation  from  the  wave  staff  was 
divided  into  the  measured  wave  height  to  obtain  the  staff  calibration  at 
that  particular  frequency.  This  calibration  was  then  normalized  by  dividing 
by  the  static  calibration.  Figure  9  summarizes  the  staff  calibration  from 
zero  to  30  Hz.  The  error  flags  become  larger  at  higher  frequencies,  be- 
cause wave  height  decreases  rapidly  with  increasing  f,  while  the  absolute 
possible  error  is  a  constant  of  the  micrometer  scale. 
2.  Processing  the  Wave  Record 

a.  Temporal  Variation  of  Mean  Squared  Wave  Height 

The  correlator,  as  wired  and  programmed  in  figure  10,  accepted 
the  wave  record,  z(t),  squared  it,  and  averaged  the  resulting  product  for 

20  seconds  on  a  continuous  basis.  The  correlator  output  was  a  voltage 

2 
representing  mean  squared  wave  height,  a  . 

n     '        i         20  sec  ?  I  9 

a   =  20"    f  z  (*)  dt   =   (z 

2 
This  voltage,  displayed  on  a  strip  chart,  revealed  that  a  itself  was 

2 
itself  a  function  of  time,  a  (t),  in  spite  of  the  20  second  integration; 

variations  on  the  order  of  20%  were  observed  over  a  ten  minute  interval 
and  were  present  for  all  wind-driven  surfaces.  The  nature  of  these  low- 
frequency  temporal  variations  was  not  known;  they  appeared  to  be  statis- 
tical, but  on  the  possibility  that  there  might  have  been  an  external 
cause,  several  attempts  were  made  to  study  this  effect. 
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Close  examination  of  the  strip  chart  revealed  that  the  fluc- 
tuations  in  c  (t)  were  of  very  low  frequency,  apparently  less  than  0.1 
Hz  (this  observation  being  consistent  with  the  20  second  time  constant 
of  the  correlator).  The  possibility  that  the  fluctuations  were  being 
introduced  in  the  electronics  was  checked  by  oscillating  the  wave  staff 

in  a  beaker  of  tank  water  as  shown  in  figure  7.  The  amplitude  of  the 

2 
oscillatory  motion  was  a  known  constant;  the  resulting  function,  a   ,  when 

displayed  on  a  strip  chart  was  also  constant,  indicating  that  the  previ- 
ously observed  variations  were  a  real,  if  unexplained,  phenomenon  of  the 
wind  driven  surfaces.  The  electrical  power  drawn  by  each  of  the  five 
centrifugal  fans  was  monitored  and  found  to  be  constant  with  time. 

As  a  second  possibility,  the  period  of  the  tank's  natural 
seiche  was  determined  by  shocking  the  system  (quickly  opening  and  closing 
the  large  drain  valve  at  the  tank's  west  end).  The  period  of  the  resulting 

water  surface  oscillation  was  3.7  seconds  —  at  least  one  order  of  magnitude 

2 
smaller  than  the  periods  of  fluctuation  observed  in  a  (t). 

The  energy  spectral  density  of  the  fluctuations  was  determined 

using  the  equipment  of  figure  11.  To  shift  the  low  frequencies  into  a 

2 
range  more  tractable  to  analysis,  the  function  a  (t)  was  tape  recorded  at 

speed  0.375  IPS,  and  replayed  for  processing  at  37.5  IPS.  It  was  found 

necessary  to  eliminate  the  DC  component  of  the  signal,  since  this  would 

produce  delta  function-like  behavior  at  f=0  on  the  spectrum  plot,  possibly 

masing  the  low  frequencies  of  interest.  DC  elimination  was  done  with  a 

highrpass,  RC  filter,  whose  parameters  are  shown  in  figure  11.  The  survi- 

2  2 

ving  signal  after  this  filter  was  a  (100t*)  -  a  (100t*)  ,  (where  t*  = 

time  in  the  playback  frame  of  reference).  The  autocorrelation  of  this 

function  was  computed  by  the  PAR-101,  then  the  cosine  transform  of  the 
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autocorrelation  was  determined  by  the  PAR-102,  which  yielded  the  energy 
spectral  density  in  accordance  with  the  Weiner-Kinchine  Theorem. 

The  resulting  plot?  shown  in  figure  11 ,  revealed  three  strong, 
low-frequency  components  with  periods  of  approximately  500,  100  and  50 
seconds.  It  is  possible  that  these  components  might  have  been  introduced 
by  the  high-pass  filter,  but  it  is  believed  that  the  fluctuations  are 
statistical,  possibly  introduced  by  non-linear  surface  wave  interaction. 
Because  of  the  behavior  of  this  variable,  it  was  decided  to  monitor  it 
continuously  during  all  scattering  experiments,  and  to  record  its  value 
at  the  instant  that  sound  scattering  data  were  taken.  This  procedure 
appeared  both  reasonable  and  consistent,  because  both  a  (t)  and  the  scat- 
tered sound  data  were  processed  through  integrators  with  identical  time 
constants.  The  mechanics  of  data  taking  are  described  in  section  III-D. 
b.  Water  Surfaces  Used  for  Scattering 

Five  water  surfaces  were  chosen  from  which  to  scatter  sound, 
covering  the  widest  possible  range  of  RMS  wave  height,  a,  and  containing 
both  swell  and  sea-swell  combination;  a  tabular  summary  of  these  surfaces 
is  given  in  figure  12.  Using  the  correlator,  Fourier  analyzer,  and 
plotter  (but  not  the  filter)  of  figure  11,  energy  spectral  density  plots 
were  made  for  the  surfaces  and  are  presented  in  figures  13  through  15. 
For  consecutive  runs  under  identical  wind  source  conditions,  slight 
shifting  of  peak  size  and  location  was  always  noted,  as  well  as  changes 
in  the  total  area  under  the  curve,  a  fact  previously  established  since 
area  under  the  curve  is  equal  to  mean  squared  wave  height,  a  (t).  Figures 
13  through  15  are  representitive  energy  spectra,  each  compiled  from  exam- 
ination of  several  spectra  for  nominally  the  "same"  statistical  surface. 
It  should  be  noted  that  nonlinear  interaction  caused  the  sea  and  swell 
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spectrum  "E"  to  be  significantly  different  from  superposition  of  the  two 
constituent  spectral  (sea  "B"  and  swell  "D"). 

A  closer  analysis  of  the  energy  spectral  density  of  wind 
driven  sea  "B"  was  made  by  plotting  the  data  on  log-log  axes.  A  continu- 
ous, four  hour  wave  record,  z(t),  was  tape  recorded,  and  replayed  at 
one-hundred-times  speed,  z(100t*),  into  a  General  Radio  Wave  Analyzer, 
model  1900-A,  using  a  bandwidth  of  3  Hz;  the  data  obtained  are  plotted  in 
figure  16.  The  curve  labelled  "noise"  represents  total  system  noise,  and 
was  obtained  by  frequency  analyzing  a  four  hour  tape  of  the  "wave  record" 
of  a  mirror  surface,  in  the  manner  described  above.  The  dependence  of 
energy  upon  frequency  at  the  high-frequency  dropoff  was  found  to  be  f"   , 
which  is  in  close  agreement  with  the  Neumann-Pierson  spectrum  for  wind 

waves  and  Phillips  theory  [Ref.  6]  which  predicts  f"  .  Scheible  [Ref.  3] 

-5  5 
obtained  f  '  for  a  slightly  smaller  sea  in  the  same  tank. 

c.  Inhomogeneity  of  the  Water  Surface 

The  wind  driven  seas  produced  were  not  statistically  homogene- 
ous. The  degree  of  non-homogeneity  in  the  RMS  wave  height  was  determined 
by  analysis  of  30  minute  wave  records  from  the  nine  surface  locations 
shown  in  figure  17.  A  wave  record  from  any  one  of  the  nine  locations  was 

processed  as  shown  in  figure  10.  The  output  of  the  correlator  was  a  30- 

2 
minute  record  of  a  (t),  which  was  tape  recorded  at  .375  IPS.  This  record 

was  then  replayed  at  3.75  IPS  into  the  analogue-tr-digital  converter  of 
the  SDS  9300  digital  computer.  The  converter  sampled  the  input  function, 
a  (10t*),  at  a  30  Hz  rate.  The  square  root  of  each  sample  was  then  ex- 
tracted, which  resulted  in  an  ensemble  of  approximately  5400  digital 
values  of  the  RMS  wave  height,  a.  The  mean  and  standard  deviation  were 
computed  for  each  of  the  nine  ensembles,  and  the  results  are  presented  in 
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figure  17.  Wind  driven  sea  "B"  was  used  throughout;  a  fetch-limited 
situation  is  indicated  by  the  significant  increase  of  wave  height  with 
fetch  distance, 

d.  Wave  Height  Distributions  for  the  Scattering  Surfaces 
As  previously  mentioned  in  section  II-A,  the  e~^  term  for 
coherent  specular  scatter  was  derived  by  Beckmann  (and  all  others)  after 
assuming  a  Gaussian  distribution  of  wave  heights  [Ref.  2].  Wave  height 
distributions  were  determined  for  the  various  scattering  surfaces  and 
are  presented  in  Appendix  B,  along  with  a  detailed  description  of  how 
the  data  were  obtained.  In  summary  of  this  appendix,  it  can  be  stated 
that  all  wind  driven  seas  —  including  the  sea-swell  combination  —  had 
height  distributions  which  resembled  the  Gaussian,  although  all  were 
skewed  to  the  right.  The  pure  swell,  however,  had  a  height  distribution 
which  resembled  that  of  a  sine  wave. 

C.  ANALYSIS  OF  THE  SCATTERED  SOUND  PRESSURE 

Continuous-wave j  rather  than  pulsed,  sound  was  used  in  the  scattering 
study;  however,,  at  each  frequency  and  for  each  transducer  position, 
pulsed  sound  was  used  to  determine  the  magnitude  of  the  unwanted  direct 
ray  between  sender  and  receiver,  as  well  as  other  non-surface  reflections, 
In  most  cases  the  ratio  of  direct-path  to  mirror-reflected  pressure  was 
about  .01,  and  in  no  case  did  it  exceed  ,03.  Unwanted  reflections  from 
the  tank  walls  were  generally  one  order  of  magnitude  smaller  than  the 
direct  path  signal '. 

1 .  Isolating  the  Coherent 8_  Received  Sound  Pressure 

It  was  assumed  that  the  sound  pressure,  p,  received  at  the  hydro- 
phone after  rough  surface  reflection  was  composed  of  two  components,  one 
coherent  with  the  transmitted  signal  and  the  other  incoherent. 
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P      ^coh   pincoh 

Use  of  the  PAR  Correlation  Function  Computer  made  possible  the  isolation 
of  the  coherent  component,  and  (in  principle)  determination  of  the  magni- 
tude of  the  incoherent  component. 

The  acoustic  pressure  corresponding  to  the  driving  voltage  was 
represented  by  Pd(t)  =  P.  cos(oot).  When  the  surface  was  in  the  mirror 
condition,  the  output  of  the  receiving  hydrophone  was  totally  coherent 
with  the  sent  signal.  After  amplification  this  received  signal  was  repre- 
sented by  BP.cos  oo(t-t'),  where  B  is  an  amplitude  scale  factor,  and  oot' 
represents  the  phase  shift  due  to  propagation  distance  and  possible  phase 
shifts  within  the  transducers.  Cross-correlation  of  sent  and  received 
signals  produced  the  function  M(x). 
B    20  sec   2 

M(t)   =   -Rjr    /      P.   COS  to(t-x)  COS  o)(t-t')  dt 

where  x  =  correlation  delay  applied  to  the  Pd(t)  input. 
The  output  of  the  correlator  was  displayed  on  an  oscilloscope,  from  which 
the  value  of  M(t)  could  be  read  as  a  function  of  the  variable  x.  The 
function  has  a  maximum  at  x  =  t'. 


M 


20  sec 


max 


=  M(t')  =  ~2TJ    f        Pd  cos  w(t-f)  dt 


0 
=  B  Pd2  cos2 (wt -oof)   =  B/pd^(t)>  =  B   Id 

Id  is  proportional  to  driving  intensity,  and  /  )  means  an  average  over 
20  seconds. 

The  water  surface  was  then  roughened  and  the  correlation  repeated, 
taking  care  not  to  change  the  driving  pressure,  p.(t).  The  scattered 
pressure  was  now  of  the  form 
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p  =  ABPdcos  (t-t1)  +  p(t). 


ncoh 


where  A  =  amplitude  reduction  factor  due  to  rough  surface.  Cross- 
correlation  of  sent  and  received  signals  then  gave  N(t),  which  was  composed 
of  two  integrals. 

20  sec 
N(x)  =  ^Jy    /    Pd2  COS  oj(t-t')  cos  oj(t-x)  dt    + 


P.   20  sec 

20  Q  ;    P(t)incoh  '  cos  ^-i)     dt 


The  second  term  vanishes  because  of  the  incoherence  of  the  two  factors 
in  the  integrand.  The  first  term  has  a  maximum  at  t  =  t'. 

ad  t-v   sec  j         ^ 
Nmax  =  N(t'>  -  W    /    Pd2cos2.(t-t')  dt 

=  AB   pd2(t)    =  AB  (Id 


Forming  the  ratio  of  N  „  to  H  ,   gives: 
3  max      max  3 


N(t) 


Ud 


B  /I, 


2 

Since  A  relates  pressures,  the  intensities  are  related  by  A  : 

I(coh)   \    I     A2 
I (mirror)  / 

where  ^coh)  =   (PcohP*coh\   and   Mirror  =  B2ld 
The  theory  of  section  II-A  predicts  that  A  =  e~9. 

2.  The  I ncoheren t  Component  of  Sea ttered  Sound  Pres sure 

In  principle  at  least9  the  ratio  of  incoherent  to  mirror  intensity 
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can  be  determined  using  the  correlator  in  both  the  auto  and  cross-corre- 
lation modes.  The  error  calculation  of  Appendix  A,  however,  shows  that 
for  small  values  of  g  the  uncertainty  involved  in  such  a  procedure  makes 
the  result  almost  meaningless;  this  situation  evolves  because  the  compu- 
tation requires  the  subtraction  of  two  nearly  equal  numbers,  each  with 
non-trivial  relative  error.  Although  not  used  for  computation  the  method 
is  outlined  briefly  below. 

With  mirror  surface,  the  received  sound  pressure  is,  as  before, 
BP.cos  ou(t-t').  Autocorrelation  of  this  signal  yields  C(t): 

B2P,2   20  sec 
C(T)  =  „i„   /    cos  o)(t-t')  cos  o)(t-t'-x)  dt 
20    0 

This  function  has  a  maximum  at  t  =  0,  equal  to 

B2P /       20  sec    ,  9 

C(r)max  =  C(0)  =   — —   /    cos^(t-t')  dt   =  B^  Id 

The  surface  is  roughened  and  the  autocorrelation  repeated;  the  received 
signal  is  ABP.cos  u(t--t')  +  P(t).jncoh  »  and  its  autocorrelation,  evaluated 

at  t  =  0  is  D(t)|t=0. 

?  2      2 
A^P  r       20  sec  7 

D(o)  =  °-~   /   cos^  (o(t-t')  dt 

20    0 

2ABP^   20  sec 


20   0 


/    P(t)incoh  cos  "(t-t')  dt 


,     20  sec  2 
"oT   0      P  ^  'incoh 


The  second  term  vanishes  upon  integration;  the  remaining  two  terms  are 

D(0)  ■   ^(jd)  +  ('incoh 
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Forming  the  ratio  D(0)/C(0)  gives  a  new  quantity  R: 


D(0) 

cloy 


A2B2 


But   B""  <  I 


\  incoh/ 
N.  mirror) 


/I(mirror) 

\  incoh  / 
/l(mirror)\ 

2 


incoh 


+  \  incoh/ 


B2  ( 


Determination  of  both  A  and  R  required  the  use  of  the  correlator; 
since  only  one  such  instrument  was  available,  the  two  quantities  could 
not  have  been  computed  simultaneously,  a  fact  which  would  have  further 
added  to  total  undertainty  in  the  final  result. 

D.  DATA  COLLECTION  PROCEDURE 

A  data  run  consisted  of  the  following  steps  performed  in  sequence: 

1.  With  mirror  surface  and  the  transducers  positioned  at  the  desired 
grazing  angle,  the  receiver  and  projector  were  physically  aligned 
for  maximum  scattered  return- 

2.  Using  pulsed  sound,  the  ratio  of  unwanted  direct-path  pressure  to 
mirror-surface-reflected-pressure  was  determined  and  recorded. 
Only  data  with  direct  path  less  than  3„0%  of  mirror  reflection 
value  were  studied, 

3.  With  mirror  surface  and  continuous-wave  sound,  the  cross- 
correlation  of  driving  voltage  with  received  signal  was  measured; 
the  value  of  the  maximum  of  this  function  was  read  from  an  oscillo- 
scope after  180  seconds  of  correlator  computation  time 
(corresponding  to  six  time  constants  of  the  integrators). 

4.  The  surface  was  then  roughened,  and  the  cross-correlation  was 
repeated.  The  surface  wave  record,  z(t),  was  tape  recorded  during 
this  correlation,  and  the  time  of  the  correlation  reading  was 
marked  by  annotating  the  tape  on  a  separate  voice  channel;  duration 
of  this  record  was  three  minutes. 
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5.  At  a  later  time  the  tape  containing  z(t)  was  replayed  at  original 
speed  and  autocorrelated  with  zero  delay  time.  The  resulting 
correlator  output  -mean  squared  wave  height  —  was  read  and 
recorded  when  directed  by  the  voice  channel,  corresponding  to  the 
time  of  previously  having  read  sound  correlation. 

The  circuitry  and  equipment  used  in  the  foregoing  steps  are  shown  in 

figures  18  through  20.  A  separate  data  sheet  was  filled  out  for  each 

data  run,  and  the  run  was  identified  by  its  most  important  parameters; 

an  example  of  run  labelling  is  given  in  figure  21. 

2 
Even  though  both  RMS  wave  height  and  A   varied  slightly  with  time 

throughout  a  data  run  (see  section  III-B-2-a)  they  were  determined  as 
nearly  simultaneously  as  possible,  each  having  been  processed  by  inte- 
grators with  the  same  time  constant  (the  same  integrators).  It  is  believed 
that  this  method  permitted  the  taking  of  meaningful  data,  the  precision  of 
which  was  not  appreciably  degraded  by  temporal  variation  of  important 
parameters. 
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IV.  ANALYSIS  OF  DATA 


A.  DEPENDENCE  OF  COHERENT,  SPECULAR  SCATTER  ON  SURFACE 
ROUGHNESS  FOR  SIMPLE  AND  COMPOUND  SEAS 

If  surface  shadowing  is  neglible,  Beckmann  and  all  others  predict 


([coh) 


=  e"^I  . ^      (Eq.  1).  There  is  no  reason,  however,  to  believe 

mirror   ^ 


that  this  function,  derived  for  a  Gaussian  distribution  of  wave  heights, 
would  govern  the  scattering  of  specular,  coherent  sound  from  pure  swell 
with  its  almost  sinusoidal  wave  record  (Appendix  B). 

To  compare  experimental  observations  with  theoretical  predictions, 

/lcoh{obs)       \      A2 

the  ratio  (  j fzzzjtz+\  /   =  — t:     was  formed  for  each  data  observa- 

\Icoh  (predict)/     e-g 

tion.  Departures  of  this  ratio  from  unity  can  be  attributed  to  one,  or 
a  combination,  of  the  following  factors: 

1.  Non-validity  of  equation  1  in  predicting  the  observed  phenomenon. 

2.  Experimental  uncertainty  (shown  in  Appendix  A  to  be  27%  for  the 
worst  case) . 

3.  Surface  shadowing. 

Figures  22  through  25  show  this  quantity  plotted  against  grazing 

anale  for  the  various  scattering  surfaces.  In  all  cases  —with  the  ex- 

/Icoh(obs)       \ 

ception  of  figure  23  -  the  ratio     (■? (nredict)/      ^s  witni'n  8^  of 

unity,  well  inside  the  worst  case  experimental  uncertainty.  Figure  24 
shows  that  coherent,  specular  scattering  from  narrow-band  swell  alone  is 
predicted  by  e~^  with  at  least  the  degree  of  precision  found  when  the 
sea  has  a  Gaussian  distribution  of  wave  heights. 

The  curve  representing  100  kHz  and  wind  driven  surface  "C"  shows  a 
marked  and  significant  departure  from  predicted  behavior  in  the  vicinity 
of  20°  grazing  angle.  Three  identical  data  runs  were  made  on  different 
days  at  20°  grazing  to  check  the  phenomenon's  reproducibility  (Figure  23); 
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all  points  were  well  above  the  predicted  value  of  unity,  and  well  beyond 
experimental  error-  This  phenomenon  was  again  manifest  during  the  study 
of  scatter  for  sound  beam  cross  versus  beam  parallel  to  sea  propagation 
direction  (see  figure  31). 

Since  this  gross  departure  from  theory  was  noted  only  in  the  vicinity 
of  e  =  20°  with  the  roughest  wind-driven  sea,  these  parameters  were  held 
constant,  and  the  frequency  was  varied  from  50  to  130  kHz  in  ten  kHz 
steps.  The  resulting  data  points,  plotted  in  figure  26,  show  the  onset 
of  the  excess  coherent  component  at  just  above  50  kHz. 

All  data  points  for  20°  grazina  angle  are  plotted  against  surface 

/Wobs5     \ 

roughness,  g,  in  figure  28.  Here,  /  j (diet)  >  appears  to  be  a 

monotonically  increasing  function  in  the  range  0  <  g  <  3.0  .  The  pheno- 
menon cannot  be  attributed  to  increasing  roughness  alone,  however,  since 
at  e  =  30°  even  higher  roughnesses  were  attained  without  this  behavior, 
as  shown  in  figure  27.  Similar  data,  taken  at  10°  grazing  angle,  are 
shown  in  figure  29;  the  phenomenon  in  question  does  not  appear,  but  it 
should  be  noted  that  for  all  points  g  <  0.4  .   Investigation  of  larger 
roughnesses  at  10°  grazing  will  require  sound  frequencies  in  the  range 
130  kHz  <  f  <  250  kHz,  a  region  in  which  the  ED0-327  transducers  were  not 
designed  to  operate   (250  kHz  is  the  high  frequency  limit  on  inputs  to 
the  PAR-101  Correlation  Function  Computer).  The  study  of  this  departure 
from  theory  is  therefore  incomplete,  and  investigation  of  larger  rough- 
nesses (through  higher  frequencies)  is  necessary, 

B.  COMPARISON  OF  CROSS-SEA  WITH  PARALLEL-SEA  SCATTERING 

The  width  of  the  anechoic  tank  limited  transducer  separation  to  52" 

(132  cm)  for  cross  sea  scatter  (figure  4).  If  the  distance  from  trans- 

2 

ducer  face  to  the  beginning  of  the  far  field  is  taken  to  be  x  =  2a  /x   , 
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where  a  =  piston  radius  =  11  cm,  then  x  =  81  cm  and  162  cm  for  50  kHz 
and  100  kHz,  respectively.,  Since  the  52"  spearation  clearly  placed  the 
ensonified  surface  in  the  near  field,  for  comparison  it  was  considered 
necessary  to  maintain  this  same  geometry  for  parallel  sea  scatter  (al- 
though much  more  space  was  available) . 

The  curves  of  figures  30  through  33,  which  compare  cross  and  parallel 
scatter,  differ  in  detail  from  the  corresponding  data  taken  at  18'  trans- 
ducer spacing,  possibly  due  to  operation  in  the  near  field.  For  this 
reason,  the  data  obtained  at  52"  separation  were  not  used  to  study  depend- 
ence of  specular,  coherent  scattering  on  roughness,  but  were  only  used  to 
infer  information  concerning  cross  versus  parallel  scattering 

/W0^     \ 

Figures  30  through  33  show  (j (diet)  >  versus  grazing  angle 

for  both  cross  and  parallel  scatter.  In  almost  eyery   case,  it  is  seen 
that  cross  scatter  gave  a  somewhat  stronger  specular  coherent  return  than 
was  observed  for  the  corresponding  parallel  case.  The  noticeable  drop 
which  occurred  at  100  kHz,  30°  grazing  angle  was  outside  the  bounds  of 
experimental  uncertainty;  since  it  was  not  present  at  the  greater  trans- 
ducer separation,  it  appears  to  be  attributable  to  either  near  field 
operation  or  to  the  reduction  in  ensonified  area  (by  a  factor  of  about  6). 

It  should  be  noted  that  the  larger  than  expected  coherent  component 
at  e  =  20c,  sea  "C"  occurs  in  figure  31,  for  both  cross  and  parallel 
scattering, 

C.  SURFACE  SHADOWING 

The  shadowing  function,  S(e),  derived  by  Wagner  [Ref.  5]  and  defined 
in  section  1 1 -B  uses  ray  acoustics,  neglecting  diffraction  as  well  as 
multiple  reflections.  The  ray  approximation  should  improve  with  increasing 
a/A  ,  becoming  precise  in  the  limit  as  a/A^<»  .  Wagner  states  that  the 
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correct  use  of  the  shadowing  function  for  specular,  coherent  intensity  is 

2 
multiplication  by  S  (e).  Thus 

lcoh_\    ,2,  .  -g    nr   /Icoh(obs)   \   -2,fl. 

— /-  S  (e)  e    >  °r  \  Irnh(predict)/-  S  (e) 

mirror/  \  cohVf      '' 


The  arguments  of  S(e)  are  e,  a,  and  L  where  L  =  surface  correlation 
distance  in  the  scattering  direction.  The  sound  wavelength,  x,  does  not 
appear,  but  use  of  the  ray  theory  implies  that  a  >>  x  . 

The  spatial  surface  correlation  function  measured  in  the  scattering 
direction  (up-down  wind)  resembles  a  damped  cosine  curve  [Ref.  7].  If 
the  correlation  distance  (distance  for  the  correlation  function  to  reach 
1/e  of  its  initial  value)  is  determined  using  the  decay  envelope,  L  14.0 
cm;  however,  if  L  is  computed  from  the  cosine- like  function  in  its  first 
quarter  period  (prior  to  the  first  crossing  of  the  zero  correlation  axis) 
L   3.0  cm  results.  The  correlation  distances  are  indicated  to  be  approxi- 
mate, since  they  were  obtained  for  a  slightly  smaller  sea  generated  by  a 
set  of  fans  which  is  no  longer  installed. 

The  transducers  used  limited  the  value  of  a/x  to  0.3  and  less. 
Although  this  obviously  does  not  satisfy  the  requirement  that  a  >>  x, 
data  obtained  were  compared  to  Wagner's  predictions  (Fig.  34)  plotted  for 
convenience  at  o/x  =  10,  Lack  of  experimental  data  in  the  region  0.5 
0,5  <<  o/x  <  10,0  precludes  extrapolation  and  leaves  unanswered  the 
question  whether  the  correlation  distance  computed  from  the  decay  envelope 
(L  -■  14.0  cm)  more  accurately  predicts  the  degree  of  shadowing  than  does 
L  =  3  0  cm 
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V,   CONCLUSIONS 

The  unpredicted  growth  of  the  coherent  component  of  scattered  inten- 
sity in  figures  23,  26,  and  28  may  be  the  same  phenomenon  described  by 
Medwin  [Ref„  1].  While  conducting  normal-incidence,  specular  scattering 
from  a  wind  driven  surface,  Medwin  observed  that  above  a  critical  fre- 
quency "the  instantaneous  scattered  intensity  is  a  nonstationary,  oscil- 
lating function  of  frequency,  varying  about  the  asymptotic  value,,"  He 
postulated  that  this  oscillation  with  frequency  was  due  to  alternate 
destructive  and  constructive  interference  of  incident  and  reflected  sound 
waves  from  statistically  prominent  facets  of  the  surface,  whose  distri- 
bution of  wave  heights  was  Gaussian,,  Medwin  suggested  that  the  dimension 
causing  the  interference  was  perhaps  the  distance  between  mean  positive 
and  mean  negative  wave  heights,  and  noted  that  the  phenomenon  was  absent 
for   a  <<  X.  He  defined  the  critical  frequency  as  that  for  which  x/4  = 
facet  separation  (the  first  interference  with  increasing  frequency). 

Non-normal  sound  beam  incidence  and  surface  anisotropy  greatly  compli- 
cate definition  of  the  facets  responsible  for  interference.  Examination 
of  figures  30  through  33  (particularly  the  100  kHz  curves  of  Fig,  31) 
indicates  little,  if  any  dependence  of  the  behavior  on  surface  anisotropy; 
thus,  if  this  is  an  interference  phenomenon,  it  suggests  that  the  gover- 
ning parameters  are  functions  of  x,  a,  and  e,  and  are  not  influenced  by 
correlation  distance. 

Comparison  of  parallel-sea  with  cross-sea  scattering  shows  that  the 
cross-sea  geometry  gives  a  coherent  component  which  is  in  general  consis- 
tently larger  than  that  obtained  for  the  corresponding  parallel -sea  case 
(Figs.  30  through  33).  This  difference  in  coherent  component  may  have 
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been  due  to  the  different  off-specular  rays  for  the  two  cases;  these 

rays,  which  are  excluded  in  theoretical  solutions,  are  inevitably  collected 

during  an  experiment  in  spite  of  the  use  of  directional  transducers. 

Specular,  coherent  scattering  from  pure  swell  (whose  wave  height 
probability  density  is  strongly  non-Gaussian)  appears  to  closely  obey 
the  e"q  prediction,  although  this  function  was  derived  for  a  Gaussian 
distribution  of  surface  wave  heights.  The  high-frequency  limit  of  the 
EDO-327  transducers  and  the  relatively  small  RMS  swell  wave  height 
(a  <_  0.15  cm)  restricted  the  attainable  values  of  g  to  0  <  g  <_  0.1  . 
For  completeness,  behavior  of  specular,  coherent  scattering  in  the  range 
0.1  <_  g  <_  5.0  should  be  studied. 

The  energy  spectral  density  for  sea-swell  combination  "E"  [Fig.  15] 
shows  that  about  20%  of  total  energy  is  carried  by  the  swell.  The  corre- 
sponding wave  height  probability  density  function,  however,  appears  as 
near  to  Gaussian  as  any  of  the  purely  wind-driven  seas  (Fig.  B-5).  Non- 
linear wave  interaction  caused  the  energy  spectrum  of  the  sea-swell 
combination  to  be  significantly  different  from  that  produced  by  super- 
position of  its  two  constituent  spectra  (see  Figs.  13,  14,  and  15). 

Wagner's  Shadowing  Correction  (SeCo  II-B),  based  on  ray  acoustics, 
tended  to  overcorrect,  due  no  doubt  to  the  partial  illumination  of 
"shadowed"  areas  caused  by  wave  diffraction.  In  figure  34  it  is  not 
possible  to  extrapolate  to  a  higher  value  of  o/a,  where  the  ray  theory 
should  be  more  valid,  in  order  to  determine  whether  the  proper  correlation 
distance  to  use  in  the  computation  of  Wagner's  correction  should  be 
derived  from  the  envelope  of  the  spatial  correlation  function. 

At  the  time  of  this  experimental  work,  transducers  were  being  fabri- 
cated to  cover  the  frequency  range  120  •_  f  ■_  240  kHz.  Investigation  of 
specular  scattering  in  this  frequency  range  is  necessary  to  clarify  many 
of  the  phenomena  observed,       gg 


APPENDIX  A  -  ANALYSIS  OF  ERROR 

Several  important  calculated  quantities  concerninq  scattered  sound 
depend  on  oscilloscope  readings  and  the  wave  staff  calibration.  It  was 
possible  to  estimate  the  degrees  of  uncertainty  in  these  constituent 
quantities,  and  their  effect  on  the  relative  coherent  intensity  and 
incoherent  intensity  compared  with  theory  is  now  considered.  The  method 
used  is  the  finite  differential  approximation  to  the  derivative. 

/lcoh{obs) 
1.  Uncertainty  in  y ^(predict) 

This  ratio,  which  compares  the  observed  to  predicted  specularly 
scattered  coherent  intensity,  is  called  I  for  brevity  in  the  analysis, 


(W0^     \    „   _A!   __    A2eq 

'  WPred1ctV      e-9 


N2  ak2a2(2  sin  e)2 
2 

where:     N  =  maximum  value  of  rough  surface  cross-correlation, 

read  from  scope 
M  -  maximum  value  of  mirror  surface  cross-correlation, 

read  from  scope. 
It  is  assumed  that  k  and  0  are  known  exactly;  sources  of  uncertainty  are 
M,  N,  and  a. 


Let:         «, 


y-    =  relative  uncertainty  in  I  due  to  N 
=  relative  uncertainty  in  I  due  to  M 


6IM 
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61 

— p-   =  relative  uncertainty  in  I  due  to  a 


The  total  relative  error  is  the  sum  of  the  contributions 


51.  4.  t  6L,  +  6IM  +  61 

total       N      M      a 


By  partial  differentiation: 


m      iL    -     2Na    -     21 
U'       3N  M2  e  N 


(2)       II    =    ^eg     -    4^ 

^'       8M  M3     e  M 


(3)      |i    =    A2  2k2  o  (2  sin  e)2  eg    =    M 

da  a 


Using  the  finite  differential  approximation  and  rearranging,  equations 
(1),  (2),  and  (3)  become: 


6IN 
(4)       -/  = 

2  si 
N 

(5)  ;m  = 

-2  6M 
M 

(6)       -f  = 

2q   6a 

a 

61  total 

2      ifl- 

6M     „  6a 
M     M  a 

u     6N    6M    „  .   6a       ..       ...        ,  ..  .  '  . 

where  -sr  »  — jgj-  »  and  —   are  the  constituent  relative  uncertainties. 

It  can  be  seen  taht  the  "worst  case"  will  occur  when  the  following  three 

conditions  are  simultaneously  satisfied: 
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a.  N  is  small 

b.  g  is  large 

c.  6M  is  negative,  while  6N  and  6a   are  positive. 

It  should  be  noted  that  the  value  of  M  is  controllable;  its  magnitude  is 
set  for  almost  full  oscilloscope  deflection  at  the  start  of  each  data  run 
by  adjusting  the  correlator  input  gains. 

For  the  "worst  case"  analysis,  typical  parameters  from  a  large  rough- 
ness data  run  were  used.  In  all  cases,  the  oscilloscope  readings  were 
judged  accurate  to  t.l  cm. 


g  =  3.0 

6N_    .1  cm 
N     2.  cm 

6M  _  -.1  cm 
M    9.0  cm 


.050 


,011 


3  a  . 40  cm 


=  2(.050  +  .011  +  .075)  =  .272  =  27.; 


I 


2.     Uncertainty  in    (jtosg|rgnt) 

■  I(mirror) 


If  the  method  of  section  III-C-2  had  been  used  to  isolate  the 
incoherent  intensity,  uncertainties  of  unacceptable  magnitude  would  have 
resulted,  as  will  no  be  shown.  Define: 

u  -  /i ( i ncoh )  \      da2       (from   section  III-C-2) 

H  r  MlmiVror)/   "   R  "  A 

2 

d   r 


T 


69 


where:    D  =  maximum  value  of  rough  surface  autocorrelation, 

read  from  oscilloscope 
C  =  maximum  value  of  mirror  surface  autocorrelation, 

read  from  oscilloscope 
By  the  method  of  finite  differentials 


Htotal     KLD"   CJ  +  ^ALM     NJ 


R  -  A2 


Using  parameters  which  occurred  during  an  actual  data  run,  (Run  25-E-100) 

g     =  -36 

R     =  .713 

A2  =  .671 


6£  6M 

C  M 


6jJ  6N_ 

D  N 


-  .1 

cm 

10 

,   cm 

.1 

cm 

5, 

cm 

=     -   .01       (negative  for  worst  case) 


,02 


<5H     _     (.713)(.03)  +  2(.671)(.03)     _     i   dc     -     145y 
"H     "     ^2 "" -     1.45     -     J^_ 

The  difficulty  arises  in  the  term  (R  -  A  );  the  two  numbers  are  nearly 
equal,  and  each  has  significant  relative  error*  The  above  calculation 
gives  a  conservative  estimate  of  the  uncertainty,  since  it  does  not  con- 
sider the  non-simultaneity  of  cross  and  autocorrelations  in  the  presence 
of  statistical  temporal  variations.  In  order  to  obtain  the  incoherent 
intensity,  it  would  have  had  to  be  calculated  or  measured  directly,  rather 
than  computed  by  subtraction. 
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APPENDIX  B  -  PROBABILITY  DENSITY  FUNCTIONS 
OF  THE  SURFACE  WAVE  HEIGHTS 

The  probability  density  of  wave  heights  was  determined  for  each  of 
the  representative  scattering  surfaces  as  follows: 

1.  With  the  wave  staff  located  at  the  center  of  the  scattering  area, 
a  ten  minute  tape  recording  was  made  of  the  wave  record,  z(t). 

2.  The  wave  record  was  replayed  at  ten  times  speed  into  the  analog- 

to-  digital  converter  of  the  CI5000/SDS9300  hy  brid  computer.  The  converter 

sampled  the  analog  signal,  z(10t*),  at  a  rate  of  333  Hz,  until  an  ensemble 

4  / 

of  10  digital  wave  heights  was  generated,  (t*  is  time  in  the  playback- 
digitize  frame) . 

3.  The  mean  and  standard  deviation  of  the  ensemble  were  computed;  the 
ensemble  mean  was  then  adjusted  to  zero  by  subtracting  the  previously 
calculated  mean  from  each  ensemble  member. 

4.  The  SDS9300  digital  computer  then  sorted  the  ensemble  values  into' 
intervals  of  width  ca,  where  a  =  ensemble  standard  deviation,  and  c  was 
a  constant  to  be  optimized. 

5.  Total  "area"  under  the  resulting  probability  density  histogram 

was  10  (the  number  of  observations).  Dividing  each  ordinate  value  by 

4 
10  then  gave  a  probability  density  whose  area  was  normalized  to  unity. 

The  appearance  (and  utility)  of  the  density  functions  obtained  was 

dependent  upon  the  interval  width  selected;  larger  widths  produced  a 

smooth  plot,  but  tended  to  suppress  all  but  the  gross  features,  while  too 

narrow  an  interval  gave  an  irregular  plot  due  to  sparsity  of  interval 

population.  Figure  B-2,  made  with  an  interval  width  =  0.05a  ,  illustrates 

the  latter  phenomenon, 


71 


For  all  wave-height  probability  densities  (except  pure  swell),  a 
standard  Gaussian  curve  normalized  in  units  of  c  was  plotted  for  compar- 
ison. These  analytical  curves  are  represented  by  solid  lines,  while 
the  observed  density  function  is  given  by  points. 

The  wave-height  probability  density  of  pure  swell  (Fig.  B-4) 
resembles  that  of  a  sine  wave,  even  though  its  frequency-domain 
structure  (Fig.  14)  is  by  no  means  a  delta  function. 
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APPENDIX  C  -  SUMMARY  TABLES  OF  DATA 


-g 

^otal 

Wobs) 

Freq. 

0 

\ 

a 

A/a 

9 

e  9 

I  . 

e"9 

mirror 

LONG 

SEPARATION: 

50 

30 

2.96 

.15 

.051 

.101 

.904 

.880 

.97 

50 

25 

2.96 

.15 

.051 

.072 

.930 

.897 

,96 

50 

20 

2.96 

.155 

.052 

.051 

.950 

.920 

.98 

50 

15 

2.96 

.14 

.047 

.024 

.976 

.973 

.98 

50 

10 

2.96 

.128 

.043 

.009 

.997 

.975 

.99 

50 

5 

2.96 

.120 

.041 

.002 

1.000 

.996 

.99 

50 

30 

2.96 

.42 

.142 

.794 

.451 

.541 

1.18 

50 

20 

2.96 

.455 

.154 

.436 

.647 

.706 

1.08 

50 

15 

2.96 

.435 

.147 

.229 

.795 

.830 

1.01 

50 

10 

2.96 

.389 

.131 

.083 

.920 

.892 

.99 

50 

5 

2.96 

.427 

.144 

.025 

.975 

.985 

1.01 

50 

30 

2.96 

.07 

.024 

.022 

.978 

.987 

.99 

50 

25 

2.96 

.063 

.021 

.013 

.987 

.993 

1.00 

50 

20 

2.96 

.05 

.017 

.005 

.995 

1.000 

.99 

50 

15 

2.96 

.031 

.010 

.001 

.999 

.996 

.99 

50 

10 

2.96 

.042 

.014 

.001 

1.000 

1.000 

1.00 

50 

5 

2.96 

.035 

.012 

^0 

1.000 

1.000 

1.00 

50 

30 

2.96 

.155 

.052 

.108 

.904 

.874 

.97 

50 

25 

2.96 

.167 

.056 

.090 

.914 

.891 

.97 

50 

20 

2.96 

.15 

.051 

,047 

.954 

.938 

.96 

50 

15 

2.96 

.15 

.051 

.027 

.973 

.959 

.98 

50 

10 

2.96 

.126 

.043 

.009 

.992 

.981 

.99 

50 

5 

2.96 

.145 

.049 

.003 

1.000 

.990 

.99 

50 

25 

2.96 

.403 

.136 

.525 

.591 

.610 

1.00 

SHORT 

SEPARATION  (Wi 

th  Wind): 

50 

12 

2.96 

.134 

.045 

.014 

.986 

.940 

.97 

50 

12 

2.96 

.455 

.154 

.161 

.852 

.818 

.95 

50 

12 

2.96 

.150 

.051 

.017 

.983 

.979 

.99 

50 

12 

2.96 

.217 

.073 

.037 

.963 

.906 

.94 

50 

30 

2.96 

.141 

,048 

.089 

.914 

.857 

,94 

50 

30 

2.96 

.429 

.145 

.825 

.438 

.543 

1,13 

50 

30 

2.96 

.086 

.029 

.033 

.967 

.935 

.94 

50 

30 

2.96 

.194 

.066 

.169 

.844 

.779 

.92 

50 

20 

2.96 

.146 

.049 

.045 

.956 

.917 

.95 

50 

20 

2.96 

.422 

.142 

.375 

.687 

.747 

1.03 

50 

20 

2.96 

.116 

.039 

.028 

.972 

.959 

.98 

50 

20 

2.96 

.216 

,073 

.098 

.907 

.870 

.96 
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Freq. 

e 

A 

0 

A/a 

9 

e"9 

total 

I    . 

Wobs> 

e"g 

mirror 

SHORT 

SEPARATION 

(Cross-wind) 

50 

30 

2.96 

.153 

.052 

.105 

.900 

.902 

.986 

50 

30 

2.96 

.467 

.158 

.980 

.375 

.554 

1.310 

50 

30 

2.96 

.120 

.041 

.065 

.937 

.972 

1.000 

50 

30 

2.96 

.225 

.076 

.228 

.796 

.879 

1.013 

50 

12 

2.96 

.172 

.058 

.023 

.977 

.981 

.992 

50 

12 

2.96 

.480 

.162 

.179 

.835 

.881 

1.017 

50 

12 

2.96 

.143 

.048 

.016 

.984 

1.000 

1.005 

50 

12 

2.96 

.251 

.085 

.049 

.951 

.968 

1.003 

50 

20 

2.96 

.153 

.052 

.049 

.951 

.974 

.983 

50 

20 

2.96 

.448 

.150 

.422 

.656 

.803 

1.151 

50 

20 

2.96 

.113 

,038 

.027 

.973 

.994 

1.012 

50 

20 

2.96 

.205 

.069 

.089 

.915 

.948 

1.007 

LONG 

SEPARATION: 

100 

25 

1.48 

.15 

.101 

.290 

.748 

.765 

.93 

100 

20 

1.48 

.14 

.095 

.166 

.848 

.832 

.95 

100 

15 

1.48 

.15 

.101 

.109 

.896 

.850 

..95 

100 

10 

1.48 

.14 

.095 

.043 

.957 

.920 

.956 

100 

5 

1.48 

.117 

.079 

.008 

.993 

.954 

.957 

100 

30 

1.48 

.40 

.270 

2.88 

.056 

.189 

.95 

100 

25 

1.48 

.405 

.273 

2.11 

.121 

.241 

.98 

100 

20 

1.48 

.42 

.280 

1.48 

.228 

.422 

1.49 

100 

20 

1.48 

.46 

.310 

1.78 

.169 

.375 

2.13 

100 

20 

1.48 

.39 

.263 

1.28 

.278 

.329 

1.23 

100 

15 

1.48 

.45 

.304 

.98 

.375 

.487 

1.36 

100 

10 

1.48 

.416 

.281 

.377 

.685 

.701 

1.02 

100 

5 

1.48 

.435 

.294 

.103 

.902 

.819 

.898 

100 

30 

1.48 

.067 

.045 

.081 

.923 

.975 

1.03 

100 

25 

1.48 

.062 

.042 

.050 

.951 

.968 

.98 

100 

20 

1.48 

.039 

.026 

.013 

.987 

.976 

.977 

100 

15 

1.48 

.042 

.028 

.009 

.998 

.987 

.988 

100 

10 

1.48 

.033 

.022 

.003 

1.000 

.988 

.98 

100 

5 

1.48 

.039 

.026 

rsj  o 

1.000 

1.000 

.991 

100 

30 

1,48 

.175 

.118 

.550 

.577 

.661 

1.07 

100 

25 

1.48 

.168 

.113 

.364 

.694 

.713 

.97 

100 

20 

1.48 

.14 

.095 

.165 

.848 

.782 

.92 

100 

15 

1.48 

.16 

.108 

.124 

.882 

.824 

.942 

100 

10 

1.48 

.143 

.097 

.045 

.955 

.917 

.96 

100 

5 

1.48 

.147 

.100 

.012 

.988 

.947 

.953 
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-g 

!total 

Wobs) 

Freq. 

e 

X 

0 

\/o 

9 

e  s 

T  . 

e"9 

mirror 

SHORT 

SEPARATION  (Wi 

th  Wind): 

100 

12 

1.48 

.162 

.109 

.082 

.921 

.856 

.923 

100 

12 

1.48 

.483 

.326 

.728 

.483 

.580 

1.233 

100 

12 

1.48 

.143 

.097 

,064 

.938 

.939 

1.022 

100 

12 

1.48 

.229 

.155 

.163 

.850 

.734 

.853 

100 

30 

1.48 

.134 

.091 

.324 

.723 

.629 

.755 

100 

30 

1.48 

.422 

.215 

3.21 

.040 

.174 

.425 

100 

30 

1.48 

.104 

.070 

.195 

.822 

.819 

.620 

100 

30 

1.48 

.214 

.144 

.826 

.438 

.516 

.600 

100 

20 

1.48 

.159 

.107 

.213 

.808 

.722 

.900 

100 

20 

1.48 

.441 

.298 

1.64 

.194 

.321 

1.58 

100 

20 

1.48 

.106 

.072 

.095 

.910 

.885 

.967 

100 

20 

1.48 

.233 

.157 

.457 

.633 

.662 

.905 

SHORT 

SEPARATION  (Cross-wind) 

100 

12 

1.48 

.166 

.112 

.086 

.918 

.911 

.982 

100 

12 

1.48 

.498 

.336 

.774 

.461 

.631 

1.27 

100 

12 

1.48 

.170 

.115 

.090 

.914 

.974 

1.015 

100 

12 

1,48 

.249 

.168 

.193 

.823 

.897 

1.029 

100 

20 

1.48 

,161 

.109 

.218 

.804 

.772 

.965 

100 

20 

1.48 

.449 

.304 

1.705 

.181 

.361 

1.564 

100 

20 

1.48 

.081 

.055 

.055 

.946 

.958 

.981 

100 

20 

1.48 

.228 

.154 

.438 

.645 

.776 

1.010 

100 

30 

1.48 

.140 

.095 

.354 

.701 

.622 

.756 

100 

30 

1.48 

.429 

.290 

3.32 

.036 

.312 

.764 

100 

30 

1.48 

.066 

,045 

.079 

.924 

.958 

.870 

100 

30 

1.48 

.169 

.114 

.515 

.597 

.654 

.708 

81 
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